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I.  SUMMARY  OF  WORK  AND  DESCRIPTION  OF  WORK  ACCOMPLISHED 


The  goal  of  this  project  has  been  the  characterization  of  non-bonding  interactions 
involving  primarily  the  aluminum  atom,  and  secondarily  the  boron,  oxygen,  and  carbon 
atoms,  in  their  ground  and  lower  excited  electronic  states,  with  the  hydrogen  molecule, 
and  also  rare  gases  and  other  molecular  species.  These  interaction  potentials  are  critical 
in  the  theoretical  modeling  of  cryogenic  media,  in  particular  solid  hydrogen,  doped  with 
isolated  atoms  and  in  the  interpretation  of  their  electronic  spectra,  which  are  the  primary 
diagnostic  tool  in  the  experimental  investigation  of  such  systems.  An  understanding  of 
these  interactions  and  how  they  affect  the  structure  and  stability  of  the  doped  solid  and 
the  electronic  spectrum  of  the  dopant  chromophore  is  greatly  facilitated  by  a 
collaborative  theoretical  and  experimental  approach.  The  long-standing  and  successful 
collaboration  of  Dagdigian's  and  Alexander's  groups  illustrates  the  utility  of  a  joint 
experimental-theoretical  research  program. 

With  the  support  of  this  grant,  we  have  succeeded  in  observing  the  electronic 
fluorescence  excitation  spectra  of  a  number  of  important  target  complexes,  including 
AI-H2,  AI-N2,  AI-CH4,  AlNe,  CAr,  and  Al'^-Ar.  In  addition,  we  have  observed  higher 
excited  electronic  states  of  the  previously  observed  B-H2  and  BAr2  complexes.  These 
observations  have  allowed  us  to  make  incisive  inferences  into  the  nature  of  the  ground 
and  excited  state  potential  energy  surfaces,  in  conjunction  with  the  collaborative 
theoretical  effort.  In  particular,  we  have  determined  the  binary  interaction  potentials 
governing  the  ground  electronic  states  of  the  AI-H2,  AlNe,  and  O-H2  system  and  have 
computed  improved  potential  energy  surfaces  for  B-H2.  The  potential  energy  surfaces  of 
these  and  other  systems  have  been  used  to  predict  the  structure  and  electronic  spectrum  of 
larger  complexes.  In  particular,  the  spectrum  of  the  BAr2  complex  for  both  the  3s  <— 
2p  and  2s2p'^  2p  transitions  has  been  computed  and  found  to  agree  well  with 
our  experimental  observations. 

The  equilibrium  properties  of  a  boron  impurity  trapped  in  solid  para-hydrogen 
have  been  investigated  with  path-integral  molecular  dynamics,  using  our  derived  B-H2 
potential  energy  surfaces.  In  groundbreaking  work,  the  effect  of  the  singly-filled  nature 
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of  the  boron  2p  orbital  was  assessed  be  comparison  with  similar  simulations  in  which  an 
orientation-averaged  B-H2  potential  energy  surface  was  used.  To  help  in  the 
interpretation  of  electronic  absorption  spectra  taken  by  M.  Fajardo  and  co-workers  at  the 
Air  Force  Research  Laboratory  (Edwards  AFB,  CA),  3^  <r-  2p  absorption  spectra  of  B- 
doped  solid  para-hydrogen  were  computed  by  a  semiclassical  Franck-Condon  technique. 
These  calculations  provided  a  reasonable  rationale  for  a  spectral  assignment  and 
estimated  band  intensity. 

The  following  sections  summarize  in  more  detail  these  research  accomplishments 
and  present  some  conclusions  relevant  to  experiments  probing  B  and  A1  atoms  in 
cryogenic  solids. 

A.  Spectroscopic  Characterization  of  Non-bonding  Interactions  of  Group  niA 
Atoms 

1.  The  AlNe  Complex.  Prior  to  our  work,  only  the  heavier  Al-rare  gas 
complexes,  and  not  AlNe,  had  been  observed.  In  order  to  motivate  study  of  the  AI-H2 
complex,  we  first  sought  to  generate  the  weakly  bound  AlNe  complex.  In  this  work, 
transitions  out  of  both  spin-orbit  manifolds  of  AlNe(X^n)  to  the  C^A  and  D^II  states, 
which  correlate  with  the  Al(3^0  +  Ne  asymptote,  and  the  state,  which  correlates 
with  A1(55)  +  Ne,  were  observed  by  laser  fluorescence  excitation.  Spectroscopic 
parameters  and  dissociation  energies  of  the  ground  and  excited  electronic  states  were 
determined.  As  an  example  of  the  importance  of  considering  the  importance  of  electron 
spin  in  determining  the  strength  of  non-bonding  interactions,  we  found  that  the 
dissociation  energies  Dq  for  the  X^fl  spin-orbit  manifolds  were  significantly  different 
[14.4  ±  0.3  and  32.3  ±  0.3  cm'l  for  and  2113/2,  respectively]. 

The  potential  energy  curves  of  the  AlNe(X2n,  states,  which  correlate  with 
the  ground  Al(3p)  +  Ne  asymptote,  were  also  calculated  by  ab  initio  methods.^  When 
consideration  of  the  electrostatic  +  spin-orbit  interactions  were  taken  into  account,  the 
calculated  Dq  values  for  the  X2ni/2  and  X2n3/2  manifolds  were  in  excellent  agreement 
with  the  experimental  values. 
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2.  The  AI-H2  Complex.  We  have  employed  laser  fluorescence  excitation  to 
characterize  the  electronic  spectra  of  AI-H2  and  other  weakly-bound  Al-molecule 
complexes.  In  other  laboratories,  attempts  to  observe  the  AI-H2  complex  by  resonance- 
enhanced  multiphoton  ionization  (REMPI)  detection  schemes  failed  because  of  the  rapid 
pre-dissociation  of  the  electronically  excited  state.  In  recent  work,  several  laboratories 
have  applied  zero-kinetic  energy  electron  spectroscopy  to  obtain  vibrational  information 
on  the  ground  electronic  states  of  the  several  neutral  and  ionic  Al-molecule 
complexes.^’^  However,  this  technique  has  not  been  applied  to  the  very  weakly  bound 
AI-H2  complex;  in  any  case,  this  approach  does  not  provide  structural  information  on  the 
excited  neutral  electronic  states. 

In  our  laboratory,  we  have  employed  a  unique  detection  scheme  suitable  for  laser 
fluorescence  detection.  In  an  electronically  excited  atom-molecule  complex  and  in 
contrast  to  atom-rare  gas  binary  complexes,  we  can  anticipate  that  there  will  be  a  high 
probability  that  the  excited  state  will  decay  by  predissociation  through  coupling  with 
lower  PES's.  We  have  succeeded  in  detecting  A1  atom-molecule  complexes  by 
excitation  to  excited  electronic  states  higher  than  the  lowest  [Al(45)-M  here].  In  this 
case,  the  transition  can  be  detected  even  if  predissociation  occurs  so  long  as  the  final  state 
is  an  emitting  excited  state. 

With  this  excitation  and  detection  scheme,  we  have  observed  transitions  from  the 
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ground  vibronic  level  of  the  AI-H2  complex  to  bound  vibrational  levels  of  Al(55')-H2. 
These  transitions  were  observed  both  by  detection  of  emission  from  lower  A1  atomic 
states  (45, 3d)  and  chemiluminescence  of  AlH(Aln)  product  by  reaction  within  the 
complex.  Resonance  fluorescence  from  excited  Al(55)-H2  levels  was  very  weak.  From 
the  relative  intensities  for  formation  of  the  various  emitting  products  and  the  estimated 
decay  rate  of  the  5s  3p  transition,  the  excited-state  lifetime  was  estimated  to  be  <  35 

ps. 

Accompanying  this  experimental  study  of  the  AI-H2  complex,  we  have  carried 
out  extensive  theoretical  studies,  which  will  be  described  in  Sec.  IIIB  below.  As 
discussed  there,  the  AI-H2  complex  can  exist  in  two  nuclear  spin  modifications.  Because 
of  its  larger  binding  energy,  we  believe  that  we  are  observing  the  AI-0H2  isotopomer. 
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For  some  bands,  partial  rotational  resolution  was  achieved,  despite  the  observable 
Lorentzian  broadening.  Analysis  of  this  structure  would  confirm  the  assignment  of  the 
isotopomer  and  also  provide  an  estimate  of  the  AI-H2  separation  and  equilibrium 
geometry.  As  discussed  in  detail  in  Sec.  IIIB,  the  equilibrium  geometry  for  the  lowest 
Al(3p)-H2  PES  is  r-shaped.  Unfortunately,  we  have  not  yet  worked  out  a  definitive 
rotational  analysis,  and  this  will  be  a  subject  for  future  work. 

The  AlH  A-^X  chemiluminescence  from  reactive  decay  of  Al(55)-H2  was 
spectrally  resolved."^  The  Alfl  state  is  orbitally  degenerate;  each  rotational  level  splits 
into  a  closely  spaced  pair  of  A-doublets  of  A'  and  A"  reflection  symmetry  in  the  plane  of 
rotation.  The  A-doublets  can  be  spectrally  distinguished  since  they  separately  emit  in  the 
PIR  and  Q  branches,  respectively.  We  find  that  the  excited  A1H(A)  product  is  formed 
predominantly  in  A'  and  not  A"  A-doublet  levels.  This  A-doublet  propensity  implies  that 
the  reflection  symmetry  is  conserved  in  the  nonadiabatic  Al(55)-H2  — >  AlH(A^n)  +  H 
reaction. 

In  more  recent  work,  we  have  extended  our  study  of  the  electronic  spectrum  of 
AI-H2  to  the  transitions  to  the  3d,  4p,  and  Ad  levels.^  As  for  the  5s  <—  3p  transition,  these 
transitions  were  detected  by  monitored  emission  from  lower  A1  atomic  states  formed  by 
predissociation  of  the  excited  AI-H2  complex.  The  transition  to  the  4p  state  is  forbidden 
in  the  isolated  A1  atom  and  is  allowed  by  the  distortion  of  the  A1  atomic  wave  function  in 
the  complex.  The  Lorentzian  widths  of  the  bands  vary  widely,  in  some  cases  allowing 
resolution  of  the  rotational  lines.  The  mechanism  of  predissociation  is  investigated 
theoretically  through  ab  initio  calculation  of  C2v  cuts  of  the  excited  PES's.  We  conclude 
that  predissociation  occurs  through  coupling  with  the  repulsive  Al(4s)-H2  PES. 

3.  The  Al-N 2  Complex.  We  have  also  detected  the  AI-N2  complex  with 
vibrational  and  partial  rotational  resolution  by  excitation  to  electronic  states  correlating 
with  the  A1(55,  4d)  levels  and  detection  of  emission  from  lower  A1  atomic  levels.^  In 
previous  work  on  this  complex.  Brock  and  Duncan  investigated  this  complex  by 
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photoionization  spectroscopy.  Quantum  chemical  calculations  by  Chaban  and  Gordon 
predict  that  the  equilibrium  geometry  of  the  Al(3p)-N2  complex  to  be  linear  and  provide 
an  estimate  of  the  ground-state  dissociation  energy.  The  different  geometries  of  the 
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A1-H2  and  AI-N2  complex  are  consistent  with  the  opposite  signs  of  the  electric 
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quadmpole  moments  of  the  molecules. 

Partial  rotational  resolution  was  achieved  in  vibrational  bands  of  the  55  <—  3p 
transition,  Consistent  with  the  linear  equilibrium  geometry,  the  profiles  were  well 
reproduced  by  simulations  of  transition  and  a  Voigt  profile  consisting  of  a 

Lorentzian  component  due  to  nonradiative  decay  and  a  Gaussian  laser  line  shape  profile. 
The  fitted  Lorentzian  widths  corresponded  to  a  decay  rate  of  7  x  10^*^  s~l.  From  the 
derived  excited-state  rotational  constant,  the  5s  <-  3p  transition  is  estimated  to  involve  a 
contraction  of  the  Al-N  bond  length  by  0.6  A. 

4.  The  AI-CH4  Complex.  We  have  also  recorded  fluorescence  excitation  action 
spectrum  for  the  5s  <—  3p  and  Ad  4-  3p  transitions  in  the  AI-CH4/CD4  complex.  As  in 
the  other  A1  atom-molecule  complexes  studied,  a  progression  in  the  van  der  Waals 
stretch  mode  was  observed.  The  widths  of  the  bands  are  much  greater  than  for  AI-H2 
and  AI-N2  and  are  indicative  of  an  increased  nonradiative  decay  rate  and/or  more 
complicated  rotational  structure.  The  widths  of  the  corresponding  AI-CD4  bands  are 
significantly  less,  suggesting  that  H-atom  motion  is  involved  in  the  nonradiative  decay. 

It  was  not  possible  to  obtain  a  definitive  vibrational  assignment  for  spectra  of  these 
complexes  and  only  lower  bounds  to  the  excited-state  binding  energies  were  derived. 

As  in  the  case  of  the  AI-H2  complex,  AlH  chemiluminescence,  as  well  as 
emission  from  lower  excited  A1  atomic  levels,  was  observed  from  excited  AI-CH4.  In 
sharp  contrast  to  the  preferential  A-doublet  formation  seen  for  reaction  within  the 
A1(55)-H2  complex,  the  AlHCA^ri)  product  from  reaction  within  Al(4f/)-CH4  was  found 
to  have  approximately  equal  populations  of  the  A'  and  A"  A-doublet  levels.  This  is 
likely  the  result  of  the  fact  that  the  plane  of  rotation  of  the  AlH  product  is  not  well 
defined  with  respect  to  planes  of  symmetry  in  the  complex. 

The  excited-state  dissociation  energies  determined  from  our  studies  of  A1 
atom-molecule  complexes  are  presented  in  Table  I  and  compared  the  corresponding  Al'*’ 


ion-molecule  binding  energies.  These  values  were  taken  from  work  from  our  and  other 

,  .  ,  .  1,4,6,10-13 

laboratones. 
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Table  I.  Binding  energies  (cm“^)  of  the  excited  Al-M  and  ionic  A1+-M  complexes. 


M 

Do(Al(5i)-M) 

Do(Al(4d)-M) 

£)o(A1+-M) 

Ne 

81  a 

At 

692 

787 

982  ±5 

Kr 

1216 

1429 

1529  ±2 

H2 

300  (+160,-10) 

472  ±52 

N2 

1218  ±10 

2705  ±165 

1924  ± 175 

CH4 

>1489 

>1824 

2120  ± 105 

5.  Chemical  Reaction  within  the  Excited  B(^D)-H2  Complex.  FoUowing  an 
earlier  study  of  the  35  <-  2p  transition  in  the  B-H2  complex/'*  we  investigated  the 
transition  to  the  excited  Islp^  '^D  valence  state.  No  fluorescence  was  observed  from  the 
excited  B(2Z))-H2  complex,  but  with  fluorescence  depletion  spectroscopy  we  were  able 
to  observe  the  transition.*^  Formation  of  BH  in  its  A^Il  and  states  through  chemical 
reaction  within  the  complex  was  detected  through  observation  of  BH  A  -  X  and  ft  -  a 
chemiluminescence.*^  All  bound  Alfl  and  all  energetically  allowed  b'^IT  levels  were 
found  to  be  significantly  populated. 

Comparing  the  fluorescence  depletion  and  chemiluminescence  action  spectra,  we 
find  that  the  excited  state  chemical  reaction  occurs  most  efficiently  in  the  lower-energy 
portion  of  the  electronic  transition.  There  are  5  PES’s  which  emanate  from  the  B(}D)  + 
H2  asymptote.  The  above  observations  were  used  to  draw  inferences  on  the  nature  of 
these  surfaces.  It  is  to  be  noted  that  the  reactive  nature  of  the  B(2D)-H2  PES’s  contrast 
shaiply  with  the  apparent  nonreactive  nature  of  the  B(35)-H2  PES.  These  results  may 
have  implications  for  the  lower  excited  AI-H2  PES,  which  correlates  with  Al(45)  +  H2; 
this  PES  may  also  be  repulsive  like  the  B(35)  +  H2  PES. 

B.  Theoretical  Characterization  of  Interactions  of  Hydrogen  with  Aluminum, 
Boron,  and  Oxygen 

In  earlier  work  supported  by  an  earlier  AFOSR  grant,  we  used  multi-reference, 
configuration-interaction  (MRCl)  methods  to  determine  the  3  adiabatic  potential  energy 
surfaces  (PES’s)  describing  the  B(2p)-H2  system.*^  Then,  following  Rebentrost  and 
Lester,**  we  used  the  matrix  elements  of  the  electronic  angular  momentum  I  to  define  the 
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transformation  into  an  approximate  diabatic  basis,  characterized  by  the  3  orthogonal 
projections  of  the  Ip  orbital  along  the  B-H2  axis.  Our  work  provided  the  benchmark  for 
later  ab  initio  investigations  by  Chaban  and  Gordon  [on  the  energetics  of  M— H2 
complexes  in  high  symmetry]  and  Bauschlicher  [on  the  fragmentation  energies  of 
B(H2)n  complexes],  both  of  which  confirmed  the  accuracy  of  our  earlier  work.  More 
recently,  our  B-H2  PES's  have  been  used  by  Niday  and  Weeks,  in  a  time-dependent  study 
of  inelastic  scattering,  and,  as  will  be  described  below,  by  Voth  and  co-workers  m  a 
path-integral  simulation  of  the  properties  of  solid  pH2  embedded  with  a  number  of  B 
atoms. 

We  fitted  the  calculated  B(2p)-H2  diabatic  PES’s  to  the  general  angular 
expansion  of  Dubemet  and  Hutson  in  terms  of  reduced  rotation  matrix  elements.  These 
expansions  were  then  used^^  in  variational  calculations  of  the  bend-stretch  wave 
functions  of  the  complex.  The  adiabatic  bender  approximation  was  used  to  gam 
additional  insight  into  the  characterization  of  the  levels  of  the  complex.  Our  research 
program  is  unique  in  that  we  have  the  capability  and  expertise  both  to  determine  accurate 
ab  initio  PES's  and,  subsequently,  to  treat  the  nuclear  dynamics  correctly  to  determine  the 
rovibronic  levels  of  the  weakly-bound  complex. 

During  the  past  three  years,  we  carried  out  similar  studies  of  the  0(^P)+H2  and 
A1(^P)+H2  systems.^^’^^  The  reaction  of  0(^P)  with  H2  has  a  barrier  >13.5  kcal/mol  in 
height.  Because  of  this  high  barrier  to  reaction  it  may  be  possible  to  stabilize  ground- 
state  oxygen  atoms  in  liquid  or  solid  H2.  For  this  reason,  doping  of  atomic  O  into  solid 
H2  has  been  suggested  as  a  possible  pathway  to  enhancing  the  specific  impulse  of 
cryogenic  propellants.  Motivated  by  this  suggestion,  Li,  Apkarian,  and  Harding  (LAH) 
carried  out^^  an  ab  initio  investigation  of  the  0(^P)H2  PES's  in  the  region  of  the  van  der 
Waals  minimum. 

In  independent  work  we  carried  out  similar  MRCI  calculations  with  a  large 
(augmented  quadruple  zeta)  basis  set.  The  transformation  of  the  two  electronic  states  of 
A"  symmetry  into  a  diabatic  basis  was  performed  equivalently  to  our  earlier  work  on 
BH2.  To  correct  for  the  incomplete  recovery  of  the  electron  correlation  energy,  we 
scaled  the  external  correlation  energy,  as  originally  proposed  by  Truhlar,  Werner  and 
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their  co-workers,^®’^^  fixing  the  scaling  factor  to  match  the  results  of  coupled-cluster 
calculations  in  linear  and  perpendicular  geometry.  Our  resulting  well  depth  was 
substantially  smaller  than  that  reported  by  LAH  (79  cm“^  vs.  108  cm“^).  In  addition, 

LAH  predict  the  deepest  well  to  occur  in  linear  geometry,  whereas  we  predict  the  deepest 
well  to  occur  in  perpendicular  geometry.  We  used  several  tests  to  confirm  to  accuracy  of 
our  calculations;  the  origin  of  the  discrepancy  with  the  earlier  work  of  LAH  is  not  yet 
clear. 

Using  our  ab  initio  PES's  in  the  subsequent  simulation  of  the  bound  states  of  the 
0(^P)-H2  complex,  we  predict  zero-point  binding  energy  of  15  cm"^  and  23  cm"^  for  the 
complexes  with  pH2  and  0H2,  respectively.  This  is  considerably  smaller  than  the 
comparable  binding  energies  for  the  5(2p)-H2  complexes  (40  cm"^  and  60  cm"^  forpH2 
and  oH2).^^  The  weak  dispersion  energy,  along  with  a  favorable  quadrupole-quadrupole 
interaction,  provides  the  binding  in  these  complexes.  The  single  p  electron  of  the  B  can 
be  easily  oriented  to  maximize  the  quadrupole-quadrupole  interaction  as  well  as  reducing 
the  Pauli  repulsion  with  the  filled  H2  ICg  orbital.  However,  the  O  atom  has  four  p 
electrons.  In  its  ground  state,  two  of  these  are  triplet  coupled,  so  that  no  matter  from 
which  direction  the  H2  moiety  approaches,  at  least  one  p  electron  will  experience  Paul 
repulsion  with  the  filled  H2  ItTg  orbital.  It  is  for  this  reason  that  the  binding  is  weaker. 

We  have  also  carried  out  similar  MRCI  calculations  for  the  AI-H2  complex.  Ab 
initio  calculations  for  this  system  had  already  been  reported  by  Gordon,  Bauschlicher  and 
their  colleagues, but  only  for  high  symmetry  geometries  (collinear  and  T-shaped). 
Our  MRCI  calculations  determined  the  PES's  over  the  complete  range  of  geometry 
necessary  to  characterize  the  binary  complex.  Calculations  were  earned  out  with 
double-,  triple-,  and  quadruple-zeta  basis  sets  and  then  extrapolated  to  the  complete  basis 
set  limit.^'^  Similar  calculations  were  carried  out  for  B-H2,  to  improve  the  accuracy  of 
ouf  earlier^^  PES's.  For  the  AI-H2  calculations,  in  order  to  obtain  a  smoothly  varying 
potential  energy  surface  over  the  full  range  of  the  configuration  space,  we  find  it 
necessary  to  carry  out,  at  each  geometry,  multi-configuration  self-consistent  field 
(MCSCF)  calculations  with  a  small  specified  subset  of  reference  configurations.  The 
orbitals  and  wavefunctions  were  then  used  as  initial  input  for  both  the  complete-active- 
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space,  self-consistent-field  (CASSCF)  and  subsequent  MRCI  calculations.  The 

calculated  well-depths  in  the  minimum  energy  geometry  (T-shaped)  were  137  cm“^  and 

220  cm“^  for  B-H2  and  AI-H2,  respectively.  The  latter  number  compares  well  with  the 

33 

best  estimate  of  227  ±  8  cm“^  of  Bauschlicher  and  co-workers. 

Determination  of  the  energies  of  the  rovibronic  levels  of  the  complex  was  done  as 
in  our  earlier  work  on  B-H2.^'*  These  calculations  are  currentiy  implemented  in  our 
Hibridon  4  code,^^  (development  of  which  was  partially  supported  by  previous  AFOSR 
grants)  for  the  investigation  of  bound-states  and  scattering.  Table  II  shows  the  calculated 
zero-point  corrected  dissociation  energies  (Dq)  for  the  B-H2  and  AI-H2  complexes. 
Although  the  AI-H2  well  depth  is  nearly  twice  as  large  as  that  of  B-H2,  the  zero-point 
corrected  dissociation  energies  are  only  ~  25%  larger.  This  is  because  the  spin-orbit 
splitting  of  A1  atom  is  so  much  larger  (112  cm“^  vs.  15.3  cm“^).^^ 


Table  n.  Zero-point  corrected  dissociation  energies  (cm  of  the  lowest  bend-stretch 
levels  of  the  B(^P)-H2  and  A1(^P)-H2  complexes.^ 

B  +  pH2  B  +  0H2  A1  +  pH2  A1  +  0H2 

39?7  59I8  504  T  743 

36±18 

33.8*1  569d  39ld 

a.  The  zero  of  energy  is  defined  by  lA(‘‘'Pif^+E2(j  =  0)  for  the  complexes  withpH2  and  by 
M(Vi/2)+H20’  =  1)  for  t^^  complexes  with  0H2. 

b.  Bauschlicher  et  al  (Ref.  ) 

d.  Binding  energies  for  H2  restricted  to  its  lowest  rotational  state  (/  =  0  for  pR2 
free  rotation  limit. 

This  past  year  we  have  begun  the  study  of  the  dependence  of  the  potential  energy 
surface  on  the  H2  bond  distance  r,  which  was  held  fixed  in  our  previous  investigations. 
The  inclusion  of  the  dependence  on  r  will  be  of  importance  in  the  interpretation  of  the 
infrared  spectroscopy,  in  a  project  in  which  we  are  collaborating  with  R.  J.  Hinde  at  the 
University  of  Tennessee. 


C.  Boron  in  Larger  Clusters  and  Matrices 

In  past  AFOSR-supported  work,  we  have  used  our  binary  B-M  interaction 
potentials  (where  M  is  either  a  noble  gas  atom  or  the  pH2  molecule)  to  simulate  the 
energetics  and  nuclear  dynamics  of  larger  clusters,  with  multiple  noble-gas  or  pH2 
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ligands.^^'^^  For  the  ground  electronic  state  of  these  complexes,  we  follow  Balling  and 

Wright‘S®  and  assume  that  the  interaction  potential  is  given  by  the  lowest  state  of  a  3  x  3 

electronic  Hamiltonian.'^^  To  model  the  stmcture  and  electronic  absorption  of  these 

complexes,  we  have  used  ’  quantum  Monte-Carlo  techniques  ’  to  determine  the 

zero-point  energy  and  an  approximation  to  the  vibrational  wavefunction  of  the  complex. 

Our  first  application  to  the  simulation  of  an  optical  excitation  spectrum  was  to  the 

3^  ^  2p  transition  in  the  BAr2  cluster.'^^  Here  the  excited  state  is  non-degenerate.  More 

recently,  Yang  and  Dagdigian'^  observed  transitions  in  the  BAr2  cluster  which  can  be 

attributed  to  the  2s2p^  '^D  <-  2s^2p  '^P  valence  transition  of  the  B-atom  chromophore. 

Size-selected  fluorescence  depletion  spectra  show  that  these  features  in  the  excitation 

spectrum  are  primarily  due  to  the  ternary  complex. 

To  model  the  electronically  excited  state  of  the  complex,  we  use  the  Balling  and 

Wright^^  approach,  modified  here  to  incorporate  the  five-fold  degeneracy  of  the  state. 

The  resulting  Hamiltonian  matrix  is  10  x  10.  Each  root  of  these  matrix  corresponds  to  a 

separate  excited  state  of  the  complex,  all  of  which  can  be  accessed  in  the  transition. 

45 

Within  this  modification,  we  invoke  the  same  Franck-Condon  approximation  which  we 

used  earlier  to  model  the  3s  <-  2p  transition.  The  underlying  diatomic  B(^D)Ar  potential 

46 

energy  curves  are  those  which  reproduce  our  previous  experimental  observations  on  the 
electronic  states  emanating  from  the  Ti(}D)  +  Ar  asymptote.  The  intensity  of  the 
transition  to  each  excited  electronic  state  was  completed  by  projecting  the  eigenvector  of 
each  excited  state  onto  the  ground  electronic  state,  and  using  standard  atomic  line- 
strength  factors.  The  simulated  spectrum  provides  a  reasonable  reproduction  of  the 
experimental  fluorescence  depletion  spectrum  and  facilitates  the  understanding  of  the 
major  features. 

In  earlier  work,  we  have  used  variational  and  diffusion  Monte-Carlo  methods  to 
simulate  the  energetics  of  small  clusters  of  B  with  Ar  and  pH2.^^’  During  the  past 
several  years,  we  have  continued  this  work,  concentrating  on  single  B  atoms  site- 
substituted  in  a  hexagonal  close-packed  lattice  of  p-H2  molecules.  The  total  interaction 
Hamiltonian  is  taken  to  be  the  pairwise  sum  of  3  x  3  B-pH2  ’  interaction  matrices  and 
pairwise  pH2-pH2  scalar  potentials.'*^  Standard  periodic  box  boundary  conditions  are 
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used.^*  We  have  written  a  fast  subroutine  for  the  calculation  of  the  potential  and  its 
gradient  for  the  interaction  of  any  number  of  B  atoms  with  any  number  of  p-^2 
molecules,  in  which  we  use"^^  a  table  lookup  and  interpolation  in  the  inverse  square  of  the 
interparticle  distances.^^ 

Our  initial  investigation  involved  variational,  followed  by  diffusion  Monte-Carlo, 
calculations.  The  variational  functions  were  based  on  a  combination  of  Jastrow  and 
pinning  functions,  similar  to  those  used  by  Cheng  and  Whaley  in  their  study  of  Li  atoms 
embedded  in  solid  H2.^^  More  recentiy,  we  have  collaborated  with  Jang  and  Voth,  in 

path-integral  simulations  ’  ’  of  B  atoms  in  solid  H2.  This  past  year  Jang  and  Voth 

22 

completed  an  important  study  based  on  just  the  spherically  symmetric  component  of 

our  B-H2  interaction  potential.  During  several  visits  to  the  University  of  Utah,  a 

graduate  student  in  our  group,  Jennifer  Krumrine,  helped  Jang  to  incorporate  oiu"  new  3  x 

27  . 

3  pairwise  B-H2  Hamiltonian,  based  on  our  most  recent  ab  initio  calculations,  into  the 
Jang-Voth  path  integral  code. 

Subsequent  simulations  of  the  optical  excitation  spectra  of  B  in  pH2  have  been 
carried  out,"^^  which  involve  also  our  ab  initio  potential  for  the  electronically  excited  B(35 
^S)H2  system.  The  electronic  degeneracy  was  treated  by  means  of  the  sum  of  pairwise  6  _ 
6  Hamiltonian  matrices,  discussed  above.  Simulations  were  performed  with  a  periodic 
box  containing  1  B  atom  and  179  pH2  molecules.  The  presence  of  a  defect  was  simulated 
by  removal  of  one  pH2  in  the  nearest-neighbor  (NN)  shell.  The  3s  <—  2p  absorption 
spectra  of  the  B  chromophore  was  predicted  by  a  semiclassical  Franck-Condon  technique 
using  path  integrals  to  sample  the  quantum  lattice  configturations.  Most  importantly,  these 
simulations  indicate  that  the  3s  <-  2p  transition  in  the  solid  is  substantially  blue  shifted 
from  its  position  in  the  free  B  atom,  because  the  excited  state  is  purely  repulsive.  These 
calculations  provide  a  reasonable  explanation  of  the  electronic  spectra  of  B  in  solid  para- 
H2  recently  reported  by  Fajardo  and  co-workers,  building  upon  earlier  work  on  B- 
doped  solidn-H2.^^’  ^  We  also  argued  that  the  second  peak  seen  in  the  experimental 
spectrum  corresponds  to  the  2s2p^  (^D)  2s^2p  (^P)  transition  in  the  solid,  which  we 

did  not  include  in  our  simulation. 
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D.  VUV  Spectroscopic  Characterization  of  the  Ar-Ar  Interaction 

In  order  to  extend  the  wavelength  range  and  hence  the  types  of  weakly  bonded 
systems  which  can  be  studied,  we  have  constructed  an  apparatus  for  the  study  of 
electronic  spectra  of  complexes  in  a  free-jet  expansion  by  VUV  LIF  detection.  Tunable 
VUV  radiation  was  generated  by  resonant  4-wave  mixing  in  Mg  vapor. 

The  first  complex  studied  with  this  apparatus  was  the  Al'^-Ar  complex.^^  The 
laser  fluorescence  excitation  spectrum  of  this  complex  was  recorded  in  the  spectral  region 
around  the  A1+  ion  3s3p  <-  3^2  iS  transition  at  167.1  nm.  A  band  system 

of  the  Al'*’— Ar  complex  was  observed.  Based  on  the  previously  reported  ground  state 
dissociation  energy  Dq"=  982  ±  5  cm“l,  Dq'  for  the  excited  state  was  determined  to 
equal  1036  ±  5  cm-1.  The  rotational  constant  of  the  state  is  close  to  those  of  neutral 
AlAr  Rydberg  states.  The  ground  and  excited  states  of  this  transition  has  electron 
occupancy  3sd^  and  3sd3pa.  The  excited  state  is  more  strongly  bound  than  the  ground 
state  despite  the  fact  that  the  repulsion  of  the  more  directional  3pG  orbital  would  be 
expected  to  be  greater  than  for  the  3s(7  orbital. 

E.  VUV  Spectroscopic  Characterization  of  the  C-Ar  Interaction 

We  have  carried  out  the  first  spectroscopic  study  of  non-bonding  interactions  of 
the  carbon  atom.^^  The  lowest  resonance  transitions  of  the  C  atom  are  in  the  VUV,  the 
2s^2p3s  <r-  2s^2p^  ^P  Rydberg  transition  at  165.7  nm  and  the  2s2p'^  2s‘^2p‘^  ^P 

valence  transition  at  156.1  nm.  Naumkin  and  McCourt  have  calculated  cb  initio 
potential  energy  curves  for  the  CAr  electronic  states  which  arise  from  the  ground  state 
C{2s^2p^  3/>)  +  At  asymptote.  The  ground  electronic  state  of  CAr  is  X^E" ,  which  has 
electron  occupancy  2s^2p7fi,  for  which  the  repulsion  is  less  than  in  the  low-lying  A^II 
state  [2scr22p;9jcr  electron  occupancy]. 

Sohlberg  and  Yarkony^^  have  carried  out  a  computational  study  of  the  binding 
and  nonradiative  decay  of  the  CAr(B^II)  state,  which  arises  from  the  first  excited  triplet 
asymptote,  C(2s'^2p3s  ^P)  +  Ar.  This  state  is  predicted  to  be  strongly  bound  [Dg  =  5100 
cm-1]  and  could  be  convenient  in  the  laser  fluorescence  detection  of  the  CAr  complex, 
but  it  might  be  predissociated  by  spin-orbit  coupling  with  a  repulsive  quintet  state.  Their 


12 


calculations  suggest  that  the  high  vibrational  levels  which  are  Franck-Condon  accessible 
from  the  CAr  ground  electronic  state  should  not  be  strongly  affected  by  this 
predissociation. 

In  our  experiments,  carbon  atoms  were  generated  in  a  supersonic  expansion  by 
193  nm  photolysis  of  several  different  precursors  [CCU,  (CH3)3CO,  CH3Br].  A  single 
weak  band,  assigned  to  the  CAr  B -X  transition,  was  observed  near  the  165.1  nm  atomic 
resonance  transition.  In  the  region  of  the  156.1  nm  transition  an  excited-state  (v',0) 
progression  was  observed.  The  vibrational  assignment  was  confirmed  by  isotopic 
substitution  of  a  CCU  photolytic  precursor,  and  the  band  system  was  assigned  as  a  - 
X^IT  transition,  with  excited  state  electron  occupancy  Isc^pT&d^.  Comparison  of  the 
partially  resolved  rotational  contours  yielded  an  excited-state  rotational  constant  Bg'  = 

0.1 1 1  cm"k  The  dissociation  energies  of  the  upper  and  lower  states  of  this  band  system 
were  determined:  Dq  =  54,8  ±  0.5  cm“^  andDo^^  =  120.4  ±  0.5  cm~^.  The  latter  is 
somewhat  larger  than  the  dissociation  energy  97.4  cm~^  obtained  from  the  potential 
energy  curves  computed  by  Naumkin  and  McCourt. 

From  the  fact  that  only  a  single  weak  band  was  observed  in  the  region  of  the 
<—  atomic  transition,  we  conclude  that  the  CAr(5^n)  state  is  significantly  affected  by 
predissociation.  We  were  also  unable  to  detect  a  dipole-allowed  transition  to  an  expected 
strongly  bound  state  correlating  with  C(2s'2p^  ^D)  +  Ar,  in  addition  to  the  weakly  bound 
state.  This  suggests  that  all  strongly  bound  states  of  CAr  may  be  predissociated  by 
repulsive  states  correlating  with  lower  atomic  states. 

F.  Observation  and  Spectroscopic  Characterization  of  AINC  and  AlCN  Isomers 

In  the  course  of  our  study  of  the  AI-N2  complex,^  we  observed  a  set  of  bands  in 
the  wavenumber  range  35  800  -  37  700  cm-l  not  assignable  to  any  previously  reported 
electronic  transitions.  These  bands  have  been  ascribed^^  to  the  AINC  molecule  from  the 
observed  change  in  the  rotational  constant  upon  isotopic  substitution  and  the 
agreement  with  the  rotational  constant  measured  by  pure  rotational  spectroscopy.^^  This 
transition  has  been  assigned  as  a  ^II  <-  XlS"*"  electronic  excitation  and  is  analogous  to  the 
Ain  <—  transition  of  the  aluminum  monohalides.  The  excited  state  is  orbitally 
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degenerate  and  is  subject  to  a  large  Renner-Teller  splitting,  as  in  the  isoelectronic  C3 
A^Iiu  electronic  state. 

We  also  observed  one  band  at  28  755  cm-^  which  we  assign  to  the  AlCN  isomer. 

It  is  interesting  that  both  the  AINC  and  the  less  stable  AlCN  isomer  are  observed  in  our 
supersonic  beam  of  photolyzed  trimethylaluminum  (T^A)  seeded  in  N2  and  a  rare  gas. 
The  observation  of  these  triatomic  species  is  surprising  since  193  nm  photodissociation  of 
N2  is  a  multiphoton  process  and  likely  not  very  efficient.  It  is  possible  that  N2  bond 
cleavage  occurs  through  reaction  with  electronically  excited  A1  atoms  formed  by 
multiphoton  absorption  of  TMA.  It  is  likely  that  supersonic  beams  containing  photolyzed 
precursors  will  be  productive  sources  of  other  novel  transient  species. 
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